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RADAR  WAVEFORMS  DERIVED  FROM  ORTHOGONAL  MATRICES 


I.  INTRODUCTION 

Generally,  modern  radars  incorporate  pulse  compression  waveforms  to  simultaneously  obtain 
sufficient  energy  on  target  for  detection  and  the  desired  range  resolution.  This  is  achieved  by  coding 
a  relatively  long  transmit  waveform  to  achieve  the  bandwidth  desired  so  that  upon  reception  and 
matched  filtering  me  pulse  is  compressed  in  time  to  a  width  that  is  approximately  equal  to  the  recipro¬ 
cal  of  the  bandwidth.  The  compressed  pulse  is  equivalent  to  the  autocorrelation  of  the  transmit  pulse 
in  the  absence  of  noise  and  Doppler  shift.  To  not  mask  weak  targets  in  the  sidelobes  of  a  stronger 
target,  the  sidelobes  of  the  compressed  pulse  should  be  small. 

Previous  pulse  compression  waveforms  are  exemplified  by  the  Frank,  Barker,  pseudo-random 
shift  register  codes,  and  the  Lewis-Kretschmer  polyphase  codes  [!].  In  the  present  report  we  describe 
new  pulse  compression  waveforms  having  low  or  zero  sidelobes  when  an  individual  waveform  or 
multiple  dissimilar  waveforms  are  processed.  These  waveforms,  also  referred  to  as  coded  sequences, 
are  related  to  certain  orthogonal  matrices  that  are  associated  with  complementary  sequences  and  also 
periodic  waveforms  having  constant  amplitude  autocorrelation  sidelobes  that  are  equal  to  zero.  We 
also  investigate  periodic  waveforms  derived  in  a  similar  manner  that  have  a  constant  sidelobe  level 
equal  to  -  1 . 

Complementary  sequences  are  multiple  code  sequences  having  autocorrelation  functions  (ACFs) 
that  sum  to  0  everywhere  except  at  the  match  point  where  the  correlation  peak  occurs.  New  orthogo¬ 
nal  processing  is  derived  whereby  the  summed  cross  correlations  of  complementary  sequences  are 
identically  equal  to  zero.  Several  codes  of  interest  are  shown  to  be  able  to  remove  a  nonmoving 
ambiguous  range  target  or  clutter  return.  Also,  subcomplementary  sequences  having  ACFs  that  sum 
to  0  beyond  a  certain  interval  are  described  and  the  theory  is  generalized. 

Periodic  waveforms  having  a  constant  sidelobe  level  equal  to  0  are  related  to  complementary 
sequences;  also,  they  are  of  interest  because  they  have  associated  aperiodic  waveforms  that  have  low 
sidelobes.  The  relationship  is  developed  between  these  periodic  waveforms  and  circular  Toeplitz 
matrices  having  inner  products  between  the  rows  equal  to  0.  Closely  related  to  these  waveforms  are 
those  that  are  associated  with  circular  Toeplitz  matrices  having  inner  products  equal  to  - 1 .  Several 
periodic  and  associated  aperiodic  codes  are  described  that  are  based  on  Number  Theory  considerations 
that  are  established  from  the  required  inner  products  of  the  circular  Toeplitz  matrices.  New  codes  are 
derived  from  the  product  of  two  periodic  codes,  and  from  a  generalization  of  the  Frank  and  the 
Lewis-Kretschmer  P4  periodic  codes;  and  also,  it  is  shown  that  certain  permutations  of  periodic  codes 
having  constant  sidelobes  of  - 1  or  0  are  also  periodic  codes  having  the  same  property.  In  addition,  a 
new  class  of  periodic  codes,  referred  to  as  reciprocal  codes,  are  briefly  described. 
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KRETSCHMER  AND  GERLACH 

n.  COMPLEMENTARY  CODES 
IIA.  Basic  Properties 

Complementary  sequences  are  coded  sequences  (complex  numbers  in  general)  having  autocorre¬ 
lation  functions  G,(fc),  which  when  time  aligned  and  added  together,  sum  to  zero  everywhere  except 
at  the  match  point.  That  is,  by  letting  the  ith  sequence  of  length  N  and  the  jth  unit  amplitude 
element  be  denoted  by  s,  • ,  the  i  th  discrete  autocorrelation  is  given  by 


Gi(k)  = 


At-* 

E  s*i,jsi,j+k 


;= i 


N+k 

E  sijs*ij-k 


i  =  i 


0<(<N-l 

-(N  -  1)  <  k  <  0 


(1) 


and  for  M  sequences,  the  sum  of  the  ACFs  is  given  by 

M  f  0  k  *  0, 

£  G'(*)  =  1  AIM,  k  =  0  (2) 

J  =  1  V. 

where  *  represents  complex  conjugation  and  k  denotes  an  integer  offset  from  the  match  point  at  zero. 
It  may  be  shown  that  at  the  match  point  the  S/N  ratio  is  maximized. 

One  of  the  early  papers  on  complementary  sequences  by  Golay  [2]  deals  with  pairs  of  binary 
coded  sequences  conceived  in  connection  with  an  optical  problem  of  multislit  spectrometry.  The  for¬ 
mulation  above  allows  for  the  more  general  case  considered  later  where  the  elements  of  the  sequences 
may  be  complex.  Golay  defines  a  set  of  complementary  series  as  a  pair  of  equally  long  binary 
sequences  (consisting  of  ones  and  zeros)  having  the  property  that  the  number  of  pairs  of  iike  elements 
with  a  given  separation  in  one  series  is  equal  to  the  number  of  pairs  of  unlike  elements  in  the  other 
series.  Golay  illustrates  this  with  the  two  sequences  (1001010001)  and  (1000000110)  having  three 
pairs  of  like  and  unlike  adjacent  elements,  four  pairs  of  like  and  unlike  elements  that  are  two  elements 
apart,  and  so  on.  A  later  comprehensive  paper  by  Tseng  and  Liu  [3]  describes  generalized  properties 
of  complementary  sets  of  binary  sequences  consisting  of  more  than  two  sequences.  Of  particular 
interest  in  Ref.  3  are  sets  of  subsequently  described  sequences  referred  to  as  mates  that  exhibit  a  cer¬ 
tain  kind  of  of  orthogonality.  In  a  later  paper,  Sivaswami  [4]  describes  multiphase  complementary 
sequences. 

We  now  discuss  basic  matrix  relationships  that  characterize  complementary  sequences.  Let  the 
rows  of  an  (M  x  N)  matrix  S  represent  rows  of  M  sequences,  each  of  length  N .  That  is. 
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The  ACFs  of  each  sequence  using  Eq.  (1)  may  be  written  as 


s  11s  1N> 

s  11S1,N-1 

+  S\2S\N > 

G2(ky. 

* 

s2\x2  A/* 

* 

s2152JV-1 

*^22^2  N  ■> 

GM(k): 

♦ 

SM\SMN’ 

* 

SM \SM.N~\ 

+  SM2SMN> 

from  which  it  is  observed  that  in  order  to  satisfy  Eq.  (2)  we  require  that  the  inner  product  of  the  first 
and  last  columns  of  Eq.  (3)  equal  0,  the  inner  product  of  the  first  and  ( n  -  l)th  columns  plus  the 
inner  product  of  the  second  and  nth  columns  equal  0.  and  so  on.  Letting  b(i)  represent  the  i th 
column  of  Eq.  (3),  we  can  state  the  required  conditions  for  complementary  sequences  as 

N-k 

Q(k)  =  £  (b(i),  b(i  +  k))  =  0.  k  =  1,  2,  •••  ,  N  -  1  (at 

i  =  i 


where  (c ,  f)  denotes  the  complex  inner  product  of  the  vectors  e  and  /  given  by 

M 

(eJ)  =  L  eJl 

.  l 


From  Esq.  (5)  it  is  noted  that  a  sufficient  condition  for  sequences  to  be  complementary  is  that  the 
columns  of  the  5  matrix  in  Eq.  (3)  be  mutually  orthogonal.  As  a  simple  example  consider  the  (2  x 
2)  Hadamard  matrix  given  by 


H{  1)  = 


1  1 
1  -1 


The  ACFs  of  the  first  and  second  rows  are  (l,  2,  1)  and  (-1,2.  - 1)  respectively,  which  sum  to  (0, 
4,  0).  Higher  order  Hadamard  matrices  having  mutually  orthogonal  columns  can  be  generated  by  the 
recurrence  relation 


H{k)  = 


H(k  -  1) 
H(k  -  1) 


H(k  -  1) 
-H(k  -  1) 


where  H(0)  =  1  and  H(  1)  is  given  by  the  example  shown  above.  The  next  higher  order  Hadamard 
matrix  H  (2)  is  given  by 


H{  2)  = 


1111 
1-1  1-1 
1  1-1-1 
1-1-1  1 


Tseng  and  Liu  [3]  describe  methods  for  generating  higher  order  complementary  sequence  matrices 
from  any  given  complementary  sequence  matrix  making  use  of  orthogonal  transformations. 
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Another  binary  complementary  matrix  can  be  generated  from  pseudorandom  shift  register  codes 
[5].  Consider  the  simple  three-element  shift  register  code  given  by  (1,  -  1,  -  1).  Forming  a  (4  x  4) 
matrix  consisting  of  the  (3  x  3)  submatrix  consisting  of  the  code  and  its  rotations,  and  a  fourth  row 
and  column  of  l’s,  we  have 


1-1-11 
-1  1-11 
-1-1  11 
1111 


which  has  mutually  orthogonal  columns  (and  rows)  and  hence  consists  of  complementary  sequences. 


Another  example  is  given  by  the  Frank  [6]  and  the  Lewis-Kretschmer  [7]  P4  code  matrices. 
The  Frank  code  matrix  has  elements  that  are  the  complex  conjugates  of  the  discrete  Fourier  transform 
(DFT)  matrix.  The  ( N  x  N)  Frank  code  matrix  F  is  given  by 


1 

o 

K 

K 

K  ■  ■ 

i 

1? 

wS 

1 

t 

f. 

< 

'  £ 
5: 
i 

1 

■  '  £ 

1 

K 

•  • 

yy{N-\)(N-\) 

where  WN  =  e  i2*,N  and  j  =  elU,ti  =  N/r-I.  For  example,  for  N  =  4,  we  have 


r  i 

Fo 

1 

1  1  I 

Ft 

1 

j  "I  -J 

f2 

1 

-1  1-1 

/L 

1 

~j  - 1  J 

(6) 


(7) 


which  is  to  be  seen  to  consist  of  mutually  orthogonal  columns  and  rows.  The  polyphase  Frank  code 
is  formed  by  chaining  or  concatenating  the  successive  rows  of  the  F  matrix.  A  Lewis-Kretschmer  P4 
code  of  length  N 2  can  also  be  written  in  terms  of  an  orthogonal  (N  x  N)  matrix.  The  P4  code 
polyphase  code  elements  are  given  by  the  phases  [1,7] 

4>(i )  —  ~(i  —  1  )2  +  7r(  i  —  I),  i  =  1 .  2,  •  •  •  ,  P  (8) 

where  P  is  the  number  of  code  elements.  The  Frank  and  P4  codes  have  several  interesting  properties 
that  arc  later  described  in  detail. 
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II-B.  Additional  Properties  of  Complementary  Sequences  of  Particular  Interest 

We  next  discuss  several  additional  properties  that  are  related  to  different  kinds  of  orthogonality. 
1.  Mates 

First  we  discuss  the  concept  of  mates  as  defined  by  Tseng  and  Liu  [3],  consisting  of  two  com¬ 
plementary  sets  of  binary  sequences. 

Definition:  A  complementary  set  of  sequences  is  said  to  be  a  mate  of  another  set  if  the  sum  of  the 
cross-correlation  functions  of  the  corresponding  sequences  in  these  two  sets  is  zero  everywhere. 

From  Ref.  3  we  obtain  the  following  theorem. 

Theorem  1:  Let  (S  S2,  ' ' '  Sp)  be  a  complementary  set  consisting  of  an  even  number  of  sequences 
where  S ,  and  S2,  S3  and  S4,  •  •  • ,  Sp  _ ,  and  Sp  are  pairs  of  sequences  of  the  same  length.  Then 


l s2,  (-SO,  S4,  (-S3).  •••,5p(-Sp_,)l 


is  one  of  its  mates  (~  denotes  a  time  reversal). 

Proof.  Let  Cx  -  denote  the  cross-correlation  function  between  x  and  the  time  reversal  of  y .  The 
theorem  is  proved  from  the  fact  that 

^xy  +  ^y(-x)  ~  ^xy  +  ^(-x)v  =  ^xy  ~  ^xy  ~  &) 

Note  that  the  proof  is  based  on  the  pairwise  sum  of  the  cross  correlations  of  two  sequences  x 
and  y  with  y  and  -x,  respectively.  Actually  this  theorem  is  more  general  and  is  not  restricted  to 
complementary  sets.  As  an  example  consider  the  two  sequences  5!  and  S2  given  by 


(1,  1,  1,  -1)  and  (1,  -1,  -1,  1). 


Their  ACFs  are 


(-1,0,  1,4,  1,0,  -1)  and  (1,  0,  1.4,  1,0,  1). 

We  next  form  the  mate  Mj  and  Af2,  which  is  given  by 

(1,  -1,  -1,  1)  and  (1  -1  -1  -1). 

The  cross  correlation  between  (S,,  A/|)  and  (S2,  M2)  is 

(-1,2,  1,  -2,  -1,0,  1)  and  (I,  -2,  -1,2,  1,0.  -1), 

whose  corresponding  elements  sum  to  zero.  This  is  potentially  useful  in  eliminating  ambiguous  range 
returns  and  for  reducing  mutual  interference  between  radars  operating  in  the  same  band. 
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2.  Cross  Correlation  Theorem 

Next,  we  derive  a  new  orthogonality  relationship  which  states  that  certain  sets  of  complementary 
sequence.!  have  cross-correlation  functions  that  sum  to  zero  by  using  all  pairwise  permutations.  Here, 
all  cross-correlation  function  permutations  are  required  in  order  that  their  sum  be  identically  equal  to 
zero.  This  differs  from  the  concept  of  mates  defined  by  Tseng  and  Liu  whereby  the  sum  of  two  cross 
correlations  at  a  time  is  equal  to  zero.  The  cross-correlation  theorem  may  be  formally  stated  in  the 
following  theorem,  which  is  proven  in  Appendix  A. 

Theorem  2:  If  the  rows  and  columns  of  an  (N  X  M)  matrix  are  orthogonal  and  all  the  columns 
except  one  sum  to  zero,  then  the  sum  of  all  cross  correlations  between  nonidentical  code  words  is 
zero. 


Examples  of  matrices  satisfying  the  theorem  are  the  circular  shift  register  matrices  uaving  an 
appended  row  and  column  of  l’s  described  previously  and  the  Frank  and  P4  matrices. 

3.  Orthogonality  of  Cross  Correlation  of  Frank  Code  Matrix  Rows 

A  still  different  orthogonality  relationship  is  presented  as  a  subset  of  Theorem  2  that  has  poten¬ 
tial  applications  in  removing  ambiguous  range  radar  returns.  Here,  the  cross  correlations  between 
codes,  represented  by  the  rows  of  a  matrix  with  a  given  separation,  sum  to  zero. 

Theorem  3:  Let  the  cross  correlations  between  rows  I  and  m  of  an  (N  X  N)  Frank  code  matrix  be 
represented  by  Clm  (n ).  Then 

N  -1 

£Clm(n)  =  0,  n  =  0,  ±1,  ±2,  ....  ±{N  -  1),  (10) 

/  =  o 


where  m  =  (I  +  r)  mod  N  and  r  -  1,  2,  .  .  .  ,  N  —  1. 

Proof:  See  Appendix  B 

As  an  example  of  removing  a  stationary  target  from  the  first  ambiguous  range  interval,  consider 
the  simplified  example  illustrated  in  Fig.  I  for  N  =  4.  First  F,,  F2,  and  F3  of  Eq.  (7)  are  transmit¬ 
ted  but  not  processed.  Then  F0,  F ,,  F2,  and  F 3  are  transmitted,  and  during  each  sweep  the  return 
signals  are  processed  with  a  filter  matched  to  the  most  recent  transmitted  code.  This  is  referred  to  as 
Channel  0  in  Fig.  1 .  By  summing  the  returns  from  the  four  sweeps  beginning  with  F0,  ambiguous 
range  targets  (or  clutter)  from  the  first,  second,  and  third  ambiguous  range  intervals  are  eliminated 
according  to  the  previous  theorem.  Moreover,  since  F0.  Fj,  F2.  and  F3  are  a  complementary  set  of 
sequences,  the  sidelobes  of  stationary  targets  in  the  unambiguous  range  intervals  sum  to  zero.  Other 
processing  channels  may  be  included  as  indicated  in  Fig.  1.  if  desired,  which  are  matched  to  the  other 
range  intervals,  and  stationary  clutter  from  the  mismatched  range  intervals  is  eliminated.  While  the 
theorem  has  been  shown  to  apply  to  the  Frank  code  matrix,  it  has  been  found  to  also  be  true  for  the 
Lewis-Kretschmer  PI,  P3.  and  P4  code  square  matrices,  and  also  the  P2  code  matrix  for  odd  square 
integers  only  (see  Ref.  1  for  a  description  of  the  PI,  P2,  and  P3  polyphase  codes). 

4.  Subcomplementary  Sequences 

Another  concept  of  interest  is  that  of  subcomplementary  sequences  [8 J .  These  consist  of  a  pair 
of  sequences  of  length  2k  r0  (k  =  1,  2,  ■  •  •)  having  ACFs  that  sum  to  zero  for  ah  shifts  equal  to  or 
greater  than  the  ACF  duration  r()  of  an  underlying  w^'erorm. 
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XMT  PULSES 


E1.J5L  .ELELJEL 


JEL 


PROCESSING  - 
INTERVAL. 


CHANNEL  0 

o 

MATCH  TO 

Fi  F2 

1 

F3 

F0  F, 

2 

r2 

r  3 

F0 

3 

F, 

F2 

F3 

Fig.  1  —  Example  of  processing  to  remove  stationarv 
clutter  and  resolve  range  ambiguity 


For  a  signal  S0  consisting  of  N  elements  and  having  an  ACF  duration  of  ±t„.  Sivaswami  |N| 
shows  that  the  2 N  sequence  5,  formed  by  concatenating  S0  with  itself.  (505,,i  and  the  2.V  sequence 
5,  formed  by  concatenating  S0  with  S0  phase  reversed  |S0<-S0)1  have  ACFs  that  when  summed 
have  zero  time  sidelobes  outside  ±r0.  Also,  the  resulting  ACF  amplitude  is  four  times  the  ACF  of 
the  sequence  S0.  The  procedure  described  above  may  be  repeated  so  that  the  new  pair  of  sequences 
Si  —  [^p^il  and  Si  =  [5j(  — 5|)]  are  four  times  as  long  as  5q  and  hence  have  four  times  the  energv. 
The  summed  ACFs  are  then  equal  to  eight  times  the  ACF  of  S0.  The  general  term  given  b\ 
Sivaswami  [8]  is  Sm  =  [Sm  .. ,Sm  _ and  5„,  =  [Sm  _  j(  —  Sm  _ | )] .  Any  coded  signal  having  desirable 
properties,  such  as  the  chirp.  Barker,  shift  register,  and  polyphase  codes  |l)  may  be  used  for  S0.  As 
an  example  let  S0  =  (1,  1.  1.  —  1 );  then 

5,  =  (1,  1.  1,  -1,  1,  1,  1.  -l) 

5,  =  (1.  1,  1,  -1,  -1.  -l,  -i.  n. 


The  ACF  of  50  is 


the  ACF  of  5  (  is 


(-1,0.  I,  4.  1.  0.  -1), 


(-1,  0,  1,4,  -1,0,  1,  S,  1,  0,  -1.4,  1.0.  -1), 

and  the  ACF  of  S  i  is 

(1,  0,  -1,  -4,  -3,  0,  3,  8,  3,  0.  -3,  -4.  -  1,  0.  1). 

Summing  the  ACFs  of  S  ]  and  S  |  we  have 

(0,  0,  0,  0,  -4,  0.  4.  16.  4.  0.  -4.  0.  0.  0,  0). 
which  is  equal  to  4  times  the  ACF  of  50. 

We  generalize  the  concept  of  subcomplementary  sequences  by  the  following  theorem. 

Theorem  4:  For  any  coded  sequence  S{)  a  subcomplementary  set  of  sequences  results  from  the 
Kronecker  product  of  Sn  and  a  matrix  consisting  of  a  set  of  complementary  sequences. 
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Proof:  See  App.ndi.x  C. 

The  matrix  in  Theorem  4  is  more  general  than  the  underlying  matrix  that  is  related  to  the  recur¬ 
sive  'elation  given  in  Ref.  8.  which  consists  of  two  rows  of  binary  elements.  The  matrix  in  Theorem 
-i  can  be  any  complementary  matrix  such  as  the  Frank  or  P4  code  matrices,  and  the  shift  register  code 
matrix  previously  discussed. 

III.  PERFECT  PERIODIC  AND  ASSOCIATED  APERIODIC  CODES 

The  investigation  of  periodic  codes  is  motivated  by  their  relationship  to  complementary  codes. 
b\  their  orthogonal  properties,  and  by  the  fact  that  many  good  aperiodic  codes  are  derived  from 
periodic  codes  having  low  sidelobe  levels.  This  is  exemplified  by  the  constant  amplitude  periodic 
Frank  code  [9],  which  has  a  constant  sidelobe  level  of  0.  referred  to  here  as  a  perfect  periodic  code 
(PPC).  We  show  in  this  section  that  the  Lew'is-Kretschmer  P4  code  and  the  generalized  Frank  and  P4 
codes  also  have  zero  periodic  sidelobes.  and  moreover  the  amplitude  of  the  aperiodic  P4  code  ACF  is 
invariant  with  rotation.  Doppler  properties  of  a  PPC  are  compared  with  the  corresponding  aperiodic 
code.  Also  a  new  PPC  is  described  that  is  formed  from  the  product  of  two  PPCs.  Relationships  are 
derived  between  PPCs  of  length  N 2  and  an  N  X  N  orthogonal  code  matrix. 

Figure  2  shows  a  periodic  waveform  consisting  of  concatenated,  constant  amplitude,  aperiodic 
codes  that  are  digitally  coded.  The  complex  values  of  a0,  a{,  •  •  •,  aN_ ,  make  up  the  code  word  of 
length  N  in  a  given  pulse  repetition  interval  (PRI).  This  code  word  is  repeated  in  succeeding  PRIs. 


AMPLITUDE 


PRI 


a0  ai  a2 


aN-1  a0  31  a2 


aN-1  a0  31  a2 


Fig.  2  —  Digitally  encoded  periodic  waveform 


A  PPC  is  defined  to  be  a  periodic  code  whose  ACr  has  zero  sidelobes  and  whose  amplitude  is 
uniform  (maximum  power  efficiency  =  1),  i.e.,  [«|!  -  ja?|  =  =  |ov_J.  An  asymptoti¬ 

cally  PPC  has  the  property  that  as  /V  —  oo  the  code's  ACF  has  zero  relative  sidelobes  and  its  power 
efficiency  is  one. 

1 1 1- A.  Properties  of  Periodic  Waveforms 

We  define  a  code  word  a  such  that  a  is  a  row  vector  of  length  ,V  and 

a  =  («„.  a  | ,  a  t.  ■  •  .  «v  i )  •  (III 
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where  an ,  n  =  0,  1 ,  2,  •  •  •  ,  N  -  1  are  the  elements  of  the  code  word.  A  periodic  code  is  one  that 
repeats  the  code  word  a  indefinitely.  Hence  if  apc  is  the  periodic  code  associated  with  a,  then 


V  =  a  a  a  "  ’ 

where  the  symbol  denotes  concatenation. 


(12) 


On  reception,  a  periodic  code  is  match  filtered  with  its  periodic  code  word.  The  output  of  the 
correlation  process  is  also  periodic  with  a  period,  N .  Hence,  the  matched  response  repeats  every  N 
unit  time  delays  as  does  the  sidelobe  response. 


We  form  an  N  x  N  circulant  matrix  A ,  based  on  each  of  the  possible  unit  time  delays  of  the 
received  code: 


a , 

a  2  ■  ■ 

'  aS- 1 

aN  —  1 

a  i  '  • 

aN~2 

aN~  2 

aN  -  1 

ao  '  ■ 

'  a,V-3 

a  2 

a  3  '  ' 

'  a0 

Note  that  the  inner  product  between  the  1st  row  and  the  m  +  1th  row,  m  -  0.  1 ,  ■  •  •  ,  N  —  1  is 
equal  to  the  output  of  the  correlation  process  at  the  mth  unit  time  delay  of  any  period.  Also  the  inner 
product  of  the  m  jth  row  and  the  w2th  row  is  equal  to  the  inner  product  of  the  w3th  and  ;»4th  row  if 
m\  ~  m2  —  my  -  /n4.  In  fact  if  we  denote  rm  to  be  the  inner  product  of  the  1st  and  m  +  1th  rows, 
then 


AA' 


ro 


r\ 

ro 
* 
r  i 


r2 

r\ 

ro 


rN  -  1 
rN -2 

r,\  -3 


(14) 


*  *  * 

rN- 1  rN- 2  rN- 3  '  '  '  r0 


where  /  denotes  the  matrix  conjugate  transpose.  Here  the  diagonal  elements  are  identical  and  equal  to 
the  matched  response,  and  the  off-diagonal  elements  are  associated  with  the  N  -  1  sidelobe 
responses. 

We  constrain  the  code  word  a  such  that 


|  a | 2  =  N  ,  (15) 

where  |  •  |  denotes  the  vector  magnitude.  We  call  a  periodic  code  “perfect"  if  all  the  code  elements 
have  equal  magnitude  and  all  of  the  sidelobe  responses  are  zero.  The  first  condition  implies  the  code 
is  100%  power  efficient  or 


a„  |  2  -  1 ,  n  =0,  1 ,  2,  •  •  •  ,  A  -  1 . 


(16) 
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Zero  sidelobe  response  implies  that 

AA‘  =  Nl,  (I7, 

where  /  is  the  N  x  N  identity  matrix. 

Equation  (17)  indicates  that  the  A  matrix  consists  of  orthogonal  rows,  and  for  a  circulant  matrix 
this  implies  that  the  columns  are  also  orthogonal  and  hence  the  A  matrix  also  consists  of  a  set  of  com¬ 
plementary  sequences.  Thus  complementary  sequences  can  be  generated  from  the  circulant  matrix  of 
the  Frank  code  and  other  codes  that  are  discussed  in  a  later  section,  and  also  from  the  Frank  code 
matrix  previously  described. 

III-B.  P4  Code 

In  Appendix  D,  it  is  shown  that  the  Lewis-Kretschmer  P4  code  has  a  periodic  ACF  with  0 
sidelobes. 

1.  Properties  of  P4  Aperiodic  Code 

We  have  shown  that  the  P4  code,  like  the  Frank  code,  is  a  PPC.  The  P4  polyphase  code  was 
originally  derived  by  sampling  the  phases  of  chirp  signal  and,  like  the  chirp  signal,  it  has  good 
Doppler  properties  and  low  sidelobes.  Figure  3  shows  a  100-element  P4  code  ACF  for  zero  Doppler. 
Figure  4  shows  the  ambiguity  surface  for  the  same  code.  The  output  of  a  receiver  matched  to  a  P4 
code  is  represented  by  Fig.  4  for  any  given  Doppler  shifted  return  signal.  Thus,  the  cut  along  the 
zero  Doppler  axis  corresponds  to  the  ACF  shown  in  Fig.  3.  A  Doppler  shifted  return  signal  causes  a 
mismatch  thereby  changing  the  filtered  signal  as  shown  by  the  ambiguity  diagram.  The  time  delay 
axis  is  normalized  to  the  uncompressed  pulsewidth,  and  the  sample  numbers  along  the  delay  axis 
correspond  to  range  cells.  The  Doppler  axis  of  the  ambiguity  function  is  given  in  terms  of  the  pro¬ 
duct  of  the  Doppler  frequency  and  the  uncompressed  pulse  duration.  A  value  of  unity  corresponds  to 
a  2 it  phase  shift  due  to  Doppler  across  the  uncompressed  pulse.  The  matched  filter  output  for  com¬ 
plex  discrete  time  samples  taken  once  per  code  element  may  be  stated  mathematically  as 


N-k 

c*(i)c(i  +  k)  e^‘  ~ 1  +*>A 

i  =  1 

0  <  k  <  N  -  1 

N+k 

£  c(i)c*(i  -  k)ejU 

i=i 

-{N  -  1)  <  k  <  0 

where 

Cj  is  the  ith  code  element, 

k  is  the  time  delay  index  relative  to  match  point,  and 
A  =  2 irfdT  (Doppler  phase  shift  in  sampling  internal  r). 


(18) 


The  ambiguity  function  is  given  by  the  squared  amplitude  of  e(k).  The  ambiguity  function  of  the  P4 
code  i.  similar  to  that  of  the  Frank  code,  which  is  shown  in  Refs.  1  and  7  to  be  derivable  from  the 
sampled  phases  of  a  step-chirp  waveform. 


2.  Invariance  of  Rotated  P4 

While  any  of  the  PPCs  remain  PPCs  if  the  underlying  aperiodic  code  is  circularly  rotated,  it  is 
not  true  in  general  that  aperiodic  codes  have  ACF  amplitudes  that  are  invariant  with  circular  rotation 
of  the  code  word.  However,  in  Appendix  E  we  prove  that  the  P4  code  has  this  invariance  property. 
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III-C.  Doppler  Properties  of  a  Periodic  Frank  Code 


Doppler  properties  of  a  periodic  Frank  code  were  investigated  to  determine  the  effect  on  the 
zero  side  lobe  level  achieved  under  zero  Doppler  conditions.  This  was  done  on  a  computer  by  consid¬ 
ering  two  concatenated  Frank  codes  as  shown  in  Fig.  5  as  the  return  signal  having  a  specified 
Doppler  frequency  fd.  The  radar  receiver  complex  video  signal  was  assumed  to  be  matched  over  one 
code  length  T  to  a  zero-Doppler  received  Frank  code,  and  the  resultant  computer  plots  were  obtained 
by  performing  the  convolution  indicated  in  Fig.  5.  Hence,  we  can  show  on  each  plot  the  resultant 
aperiodic  and  periodic  sidelobes  and  peaks.  In  Fig.  6  the  aperiodic  sidelobes  are  shown  to  the  right 
and  left  of  the  peaks,  and  the  periodic  sidelobes  are  shown  between  the  two  peaks.  It  was  found  that 
the  peak  sidelobe  levels  of  the  aperiodic  and  periodic  ACFs  are  nearly  the  same  for  Doppler  shifts  of 
0.375  and  above. 


FRANK 

FRANK 

CODE 

CODE 

0  T  2T 


MATCHED 

FILTER 


Fig-  5  —  Determination  of  aperiodic  and  periodic  sidelobes  of  a  100-element  Frank  code 


III-D.  Relationship  Between  Orthogonal  N  x  N  Matrices  and 
Perfect  Periodic  Codes  of  Length  N2 

We  address  here  the  question  of  what  the  properties  of  an  N  x  N  matrix  must  be  in  order  to  be 
able  to  concatenate  the  rows  of  this  matrix  to  form  a  code  of  length  N2,  which  when  periodically 
repeated  is  a  perfect  periodic  code.  The  sufficient  conditions  for  this  to  apply  are  stated  in  the  fol¬ 
lowing  Theorem. 

Theorem  5:  Consider  an  N  x  N  matrix  E  with  elements  on  the  unit  circle.  Let  E  have  mutually 
orthogonal  rows.  In  addition,  let  all  rotations  of  any  two  columns  of  E  be  mutually  orthogonal. 
Then  a  perfect  periodic  code  results  from  concatenating  the  rows  of  E . 

Proof.  See  Appendix  F. 

Note  that  any  N  X  N  matrix  satisfying  the  conditions  of  Theorem  5  also  has  rows  corresponding  to  a 
set  of  complementary  sequences. 

III-E.  New  Periodic  Codes 

l .  Generalized  Frank  and  P4  codes 

From  the  previous  section  we  can  generate  new  PPCs  that  we  call  generalized  Frank  codes  by 
postmultiplying  the  Frank  code  matrix  by  a  diagonal  matrix  D  consisting  of  the  elements 
(d().  </|.  d2  ‘  ’  •  ,  dN_x).  This  has  the  effect  of  multiplying  each  element  in  column  1  of  Eq.  (6)  by 
c/0,  in  column  2  by  dx,  and  so  on.  Since  the  Frank  code  matrix  meets  the  conditions  of  Theorem  5, 
and  these  conditions  are  not  affected  by  the  matrix  D,  we  see  that  the  resultant  N  x  N  matrix 
corresponds  to  a  PPC.  Also,  it  has  been  found  by  computer  simulation  that  a  P4  code  of  length  ;V: 
can  be  written  in  terms  of  an  N  x  N  matrix  satisfying  Theorem  5  and  we  can  then  generate,  in  a 
similar  manner  to  that  described  above,  a  generalized  P4  code  that  retains  the  properties  of  being  a 
PPC. 
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2.  Frank-Lewis-Kretschmer  P4  Product  Code 

Another  new  PPC  may  be  generated  by  forming  the  product,  code  element  by  code  element,  of 
odd  length  Frank  and  P4  codes.  This  is  proved  in  Appendix  G. 

3.  Permutation  Codes 

Appendix  H  shows  that  if  a(n),  n  =  0,  1,  2,  ...  N  -  1  is  a  code  having  0  or  - 1  periodic 
sidelobes,  then  aMn  ^  N,  n  =  0,  1,  ...  N  -  1,  where  M  is  an  integer  that  is  relatively  prime  to 
A,  also  has  the  same  periodic  sidelobes. 


Let  a 

= 

a  1 

°2 

«3 

a4 

Code  1 

= 

a4 

a  3 

M  =  2 

Code  2 

=  «0 

a\ 

a4 

a2 

II 

Code  3 

= 

a4 

a2 

>  M  =  4 

Codes  1,  2  and  3  have  the  same  - 1  or  0  periodic  sidelobes  as  a. 

4.  Reciprocal  Codes 

In  a  companion  report  written  by  the  authors  [10J,  a  new  class  of  codes  is  derived  that  is 
referred  to  as  reciprocal  codes.  They  are  based  on  the  fact  that  any  circulant  matrix  A  as  given  by 
Eq.  (13)  may  be  written  in  terms  of  its  eigenvectors  and  eigenvalues  as 

A  =  B  AB*,  (19) 

where  B  is  a  Butler  matrix  (discrete  Fourier  transform  matrix)  of  eigenvectors  of  A  and  A  is  a  diago¬ 
nal  matrix  of  eigenvalues  of  A .  Thus  all  circulant  matrices  have  the  same  eigenvectors  and  differ 
only  in  their  eigenvalues.  Reference  10  shows  that  as  a  consequence  of  Eq.  (19),  a  P PC  word  may 
be  multiplied  by  B  to  generate  another  PPC  word  that  consists  of  the  eigenvalues  of  the  A  matrix  of 
the  original  word. 

IV.  ASYMPTOTIC  PERFECT  PERIODIC  CODES  AND 

THEIR  APERIODIC  PROPERTIES 

An  asymptotic  perfect  periodic  code  (APPC)  is  a  periodic  code  that  becomes  perfect  as  the 
num be”  of  code  elements,  N,  in  the  periodic  code  word  approaches  infinity.  For  finite  N  either  the 
sidelobe  level  is  nonzero  and/or  the  power  efficiency  is  less  than  100%  for  these  codes.  However,  as 
A/  —  oo,  either  the  sidelobe  level  is  zero  or  approaches  zero  and/or  the  power  efficiency  is  100%  or 
approaches  100%. 

Examples  of  codes  that  are  APPCs  are  the  shift  register  code,  the  primitive  root  code,  the  qua¬ 
dratic  residue  code,  and  the  index  code.  All  of  these  codes,  except  the  index  code,  are  assumed  to 
have  a  subpulse  amplitude  of  unity,  and  have  a  sidelobe  level  (voltage)  equal  to  -1  for  all  time 
delays.  (Note  that  for  an  A/ -element  code  the  peak  or  match  point  is  assumed  to  be  N).  The  index 
code’s  first  element  is  equal  to  zero,  and  the  remaining  unit  amplitude  elements  are  polyphase. 
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Letting  e  equal  the  constant  sidelobe  level  of  a  given  APPC,  we  have  by  similar  arguments  to 
Eqs.  (11)  through  (16)  that 

N  e  6 

e  N  ■■  ■  e 

AA‘  =  .  .  (20) 

e  e  •  •  N 

The  APPCs  described  above  must  satisfy  condition  (20).  In  this  section  we  describe  these  codes 
and  also  show  the  associated  aperiodic  autocorrelation  and  ambiguity  functions.  Costas  sequences  are 
also  described  because  they  are  related  to  primitive  root  codes.  We  also  show  that  the  cross  correla¬ 
tion  between  different  primitive  root  codes  of  the  same  length  have  sidelobes  that  are  down  by 
approximately  the  time  bandwidth  product  of  the  code,  or  equivalently,  the  pulse  compression  ratio, 
which  is  determined  by  the  number  of  code  elements. 

1V-A.  Shift  Register  Codes 

These  codes,  also  known  as  pseudorandom  shift  register  codes,  when  periodically  repeated  are 
known  to  have  a  periodic  ACF  whose  sidelobes  are  a  constant  equal  to  - 1 .  These  codes  have  been 
studied  extensively  and  are  described  in  Refs.  9,  10,  and  11. 

IV-B.  Primitive  Root  Code 

The  N  code  words  of  the  primitive  root  code  [12]  are  defined  as 

an  =  Wfoi,  n  =  0,  1,2,  N  -  1 


where  N  —  1  must  be  a  prime  number  and  a  is  a  primitive  root  modulo  N  +  I  [12].  Example  plots 
of  the  aperiodic  ACF  of  the  primitive  root  code  and  its  ambiguity  function  are  given  in  Figs.  7  and  8, 
respectively,  for  N  =  100  and  a  =  2  (N  +  1  =  101  is  a  prime  number). 


Fig.  7  —  ACF  of  primitive  root  code  {p  =  101.  a  =  2) 
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Fig.  8  —  Ambiguity  function  of  the  primitive  root  code  {p  =  101,  a  =  2) 


A  different  101-element  primitive  root  code  (a  =  3)  is  shown  in  Fig.  9,  and  the  cross  correla¬ 
tion  function  (CCF)  between  the  two  primitive  root  codes  is  shown  in  Fig.  10.  These  results  are  typi¬ 
cal  for  the  cross  correlations  between  two  primitive  root  codes  of  the  same  length  or  number  of  ele¬ 
ments;  the  sidelobes  are  down  from  the  peak  by  approximately  the  pulse  compression  ratio. 


A  related  subject  that  has  received  recent  attention  in  the  literature  is  Costas  arrays  [13,14], 
This  consists  of  reordering  the  tones  of  a  step-chirp  waveform  in  order  to  eliminate  the  range-Doppler 
coupling  associated  with  the  waveform,  and  also  to  prevent  any  large  peaks  in  the  sidelobe  region  of 
the  ambiguity  function.  The  Costas  array  ordering  of  the  tones  assumes  that  no  more  than  one  tone 
will  simultaneously  correlate  with  itself  over  all  delays  and  Doppler  frequencies  in  the  pedestal  region 
of  the  ambiguity  function. 
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For  a  given  number  of  tones  A,  the  tones  are  ordered  according  to  the  sequence 

a0,  a1,  ■  ■  -  1  mod  (A  +  1), 

where  a  is  the  primitive  root  of  A  4-  1.  An  example  of  an  ACF  and  ambiguity  function  for  a  Costas 
array  is  shown  in  Figs.  1 1  and  12  for  A  =  10  and  a  =  7. 


Fig.  1 1  —  ACF  of  Costas  sequence 
(p  =  11,  a  =  7) 


IV-C.  Quadratic  Residue  Code 

For  this  binary  code,  we  introduce  the  Legendre  symbol  [12]:  {q/p).  This  symbol  is  defined 
for  all  q  that  are  not  divisible  by  p  \  it  is  equal  to  1  if  q  is  a  quadratic  residue  of  p  and  is  equal  to 
-  1  otherwise.  Note  that  q  is  a  quadratic  residue  of  p  if  the  congruence 

z2  -  q  mod  p 

has  a  solution. 
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With  these  preliminaries  the  code  is  defined  as 

a„  =  (n/N),  n  =  0,  1,  ••  • ,  N  -  1 

where  N  is  prime  number  of  the  form  4m  -  1.  Note  that  we  define  (0 /N)  — 
N  =  1 1 ,  the  code  word  a  is  given  by 

a  =  (1,  1,  -1,  1,  1,  1,  -1,  -1,  -i,  i,  -l, ). 

Plots  of  the  aperiodic  ACF  and  the  ambiguity  function  are  shown  in  Figs.  13  and 


1.  For  example  for 


14  respectively. 
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IV-D.  Index  Codes 

Let  N  be  an  odd  prime  and  a  be  a  primitive  root  of  N ,  i.e.,  {a0,  a1,  •  •  -  ,  aN~2\  mod  N  maps 
into  {I,  2,  ■■■  ,N  -  IJ.  If 

or*  =  y  mod  N , 


then  x  is  defined  as  the  index  of  y  modulo  N  to  the  base  a ,  or 

x  =  ind^y. 


In  fact  for  every  y  e{l,  2,  ■  ■  ■ ,  N  -  1)  there  is  an  associated  index  1,  •  •  •  ,N  -  2). 
For  example  for  N  =  5,  2  is  a  primitive  root  and 

2°  =1,2'=  2,  22  m  4,  23  s  3  mod  5. 

Hence 


ind2  1  =  0,  ind2  2  =  1,  ind2  3=3,  ind2  4  =  2. 

With  these  mathematical  preliminaries,  the  code  is  defined  as  (see  Appendix  I  for  derivation) 

„  ind  n  _  „  _  , ,  , 

an  =  WN_\  ,  n  =  0,  1,  2,  •  •  •  ,  N  -  1 

where  A  is  a  prime  number  and 


W, 


ind  0 


N  - 1 


=  0. 


Plots  of  the  aperiodic  ACF  of  the  index  code  and  its  ambiguity  function  are  given  in  Figs.  15  and  16 
for  N  =  101  and  a  =  3. 
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Fig.  16  —  Ambiguity  function  of  the  index  code  (101,  3) 


V.  SUMMARY 

In  this  report,  we  have  described  new  pulse  compression  waveform  coding  on  individual  and 
multiple  dissimilar  pulses  that  result  in  low  sidelobes  after  processing.  Low  sidelobes  are  important 
in  radar  applications  in  order  that  small  targets  are  not  masked  in  the  sidelobes  of  nearby  strong  tar¬ 
gets.  Also,  orthogonal  waveforms  have  been  found  that  have  cross  correlations  that  sum  to  zero 
everywhere.  This  has  potential  applications  in  resolving  ambiguous  range  targets  and  in  removing 
ambiguous  range  clutter  in  medium-  or  high-PRF  radars. 

In  particular,  we  have  extended  the  theory  of  complementary  sequences,  and  periodic  waveforms 
having  ACFs,  or  compressed  waveforms,  with  0  or  - 1  sidelobe  levels.  Complementary  sequences 
may  be  regarded  as  the  coding  modulation  that  is  applied  to  multiple  radar  waveforms.  The  comple¬ 
mentary  property  signifies  that  the  sum  of  the  individually  compressed  pulses  has  0  sidelobes. 
Periodic  coded  waveforms  having  0  or  -1  sidelobe  ACFs  are  of  interest  because  the  associated 
periodic  coded  waveforms  often  have  low  sidelobes  and,  for  the  0  sidelobe  case,  do  not  degrade  signi¬ 
ficantly  in  the  presence  of  mismatches  caused  by  Doppler-shifted  return  signals. 

By  extending  the  theory  of  complementary  sequences  and  periodic  waveforms  having  0  or  —  1 
sidelobes,  new  classes  of  waveforms  have  been  found  that  have  low  sidelobes.  Theorems  have  been 
derived  that  extend  the  existing  theory  and  that  show  new  relationships  between  complementary 
sequences  and  periodic  waveforms  having  0  sidelobes  Also,  multiple  orthogonal  coded  waveforms 
have  been  found  that  have  cross  correlation  functions  that  sum  to  zero. 
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Appendix  A 

CROSS  CORRELATION  THEOREM 


We  define  N  code  words  with  M  elements  each  as 

a „  =  (a ii  -  fli:,  m]M),  n  =  1,2 ,  •  ,  A. 


We  arrange  these  code  words  to  be  the  rows  of  the  following  N  X  M  matrix: 


a  11  a  12  '  a \M 
a2\  a22  a2  M 


A  = 


a,\\  aN2  '  aNM 


(Al) 


(A2) 


We  show  that 

Cross  Correlation  Theorem:  If  the  rows  a  A  columns  of  A  are  orthogonal  and  all  columns  except  one 
sum  to  zero,  then  the  sum  at  all  cross  correlations  between  nonidentical  code  words  is  zero. 

Proof.  A  correlation  vector  of  length  2M  -  1  between  code  words  is  defined  as 


=  am  *  a„ 


m  ,  n  =  1.2, 


,N 


(A3) 


where  *  denotes  linear  convolution  and  -  denotes  the  time-reversed  complex  conjugate  of  a„ .  The 
center  element  of  cm„  is  called  the  match  point.  This  point  can  be  computed  by  taking  the  inner  pro¬ 
ducts  between  the  rows  of  A.  However,  since  the  rows  are  orthogonal,  the  contributions  from  the 
cross  correlations  (m  *  n )  to  the  match  point  is  zero.  Hence 


Match  Point 


£ 


m  .n  -  1 
m  ±  n 


o. 


(A4) 


If  the  columns  are  orthogonal,  it  is  known  that  the  codes  are  complementary  or  that 


f  N 

Sidelobes  £  c„„ 

l"  - 1 


=  0. 


(A5) 
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Appendix  B 

PROOF  OF  THEOREM  3 


Theorem  3:  Let  the  cross  correlations  between  rows  l  and  m  of  an  (N  xN)  Frank  code  matrix  be 
represented  by  Ctm  (n ).  Then 


V  Clm(n)  =  0,  n  =  0,  ±1,  ±2 . ±(N  - 1) 

1=0 

where  m  =  (l  +  r)  mod  N  and  r  —  1,2 ,  .  . .  ,  (N  —  1). 

Proof.  The  cross  correlation  function  of  the  two  complex  sequences  x(k)  and  y(k),  each  of  length  N, 
is  given  by 


Clm (n )  =  £  x(k)y*(k  -  n)  , 

k  =0 
N-n-  1 

C,m(n)  =  x(k  +  n)y*(k)  , 

k=  0 

Let  x(k)  and  y(k),  k  =0,  1,  .. .  ,  N  -  1,  denote  the  / th  and  with  rows  of  a  Frank  code  matrix, 
where 


for  n  <  0  , 

for  n  >  0  . 


(Bl) 


kl 

J 2*T7 

x(k)  =  e  'v  , 
y(k)  =  e 


(B2) 


The  cross  correlation  function  between  rows  /  and  m,  Clm(n ),  where  /  =0,  1 , 
m  =  (l  +  r )mod N  is  given  by 


N  -  1  and 


nm  n  kit  -  m) 

jU~—  Ntn- I  j  2ir - - - 

C,m  (n)  =  e  N  2  e  N 

k  =0 

,2  N-n-l 

Clm(n)  =  e  N  £  e 

k  =o 


for 

for 


n  <  0  , 
n  >  0  . 


(B3) 
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For  m  =  (/  +  r )  mod  N . 


/2,n('+r)  N+n- 1  ->2,*r 

C/M(n)  =  *  "  L  e  N, 

for 

n  <  0  , 

k  =0 

y2TA/  N~n~l 

Clm(n)  =  e  N  £  *  "  , 

for 

n  >  0  . 

(B4) 

=0 

From  Eq.  (B4)  it  is  elementary  to  show  that 

N-l 

S  C/m  («)  =  0  ,  rt  =  0,  ±  1,  ±2,  . . 

•  ,  ±(N  - 

1)  Q.E.D. 

(B5) 

1=0 

i 
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Appendix  C 

PROOF  OF  THEOREM  4 

Theorem  4:  For  any  coded  sequence  S0,  a  subcomplementary  set  of  sequences  results  from  the 
Kronecker  product  of  S0  and  a  matrix  consisting  of  a  set  of  complementary  sequences. 


Proof.  Consider  the  Kronecker  product  of  any  coded  sequence  S0  and  a  complementary  set  of 
sequences  that  comprise  the  matrix  B  given  by 


11 

b  12  • 

"  b\N 

21 

b  22 

'■  b2N 

B  = 


(Cl) 


b\f  i  bM  2  bMN 


The  Kronecker  product  of  S0  and  B  results  in  the  i  th  row  being 

(bnS0  bj2S0  bjjSo  •••  biNS0).  (C2) 


Letting  X(t)  denote  the  autocorrelation  function  of  the  L  -element  code  50,  the  ACF  of  the  sequence 
given  by  Eq.  (C2),  which  we  denote  by  G,(f ),  is  easily  shown  to  be 


G,(t)  =  (bn  b*N)  X  (t)  (C3) 

+  (bj i bj  N _ i  +  bj2biN)  X  (t  —  L) 

+  (bi i b’,N -2  +  bj2 b*jj _ ]  +  bi3b*N)X  (t  -2 L) 

+  (b, ,  b’{  +  bi2b*2  +  •  ■  •  biN  b*N)  X  ((t  -  {N  -  1  )L] 

+  (£>;,  bN)  X  [t  -  2(N  -  1  )L}. 

From  Eq.  (C3)  it  is  observed  that  the  coefficients  of  the  X(-)  terms  are  exactly  the  same  as  the  ACF 
of  the  (th  row  of  the  complementary  matrix  B.  Hence,  it  follows  from  the  property  of  complemen¬ 
tary  sequences  that 

M 

^  G,(t)  =  NM  X[t  -  (N  -  1)L).  (C4) 
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DERIVATION  OF  THE  P4/CHU  CODE 


In  this  appendix,  we  derive  the  Chu  code  [Dl]  by  using  the  concept  of  integer  functions.  Chu 
stated  this  code  without  a  derivation,  then  showed  that  it  had  periodic  sidelobes  equal  to  0.  The 
Lewis-Kretschmer  P4  code  (8)  and  its  permutation  set  (see  III  E3)  are  the  same  as  the  Chu  code 
except  for  a  linear  phase  shift  that  does  not  affect  the  magnitude  of  the  periodic  or  aperiodic  auto¬ 
correlation  function.  Consider  the  N  x  N  circulant  matrix,  A  ,  associated  with  a  periodic  code.  Let 
the  first-row  elements  have  the  following  form: 

an  =  Wl{n\  n  =  0,  1,  2,  •  •  ■,  N  -  1  (Dl) 

where  /(■)  is  some  real  function  of  the  integer  n  and 

_  •  2* 

WN  =  e  J  N  .  (D2) 

Note  that  =  1.  We  assume  the  following  constraint  equation  (congruence)  on/(  ): 

f(n )  s  f(n  +  N)  mod  N  .  (D3) 

This  congruence  indicates  the  periodicity  of  /  and  implies  that  f (n)  -  / (n  +  N)  is  an  integer  that  is 
divisible  by  N . 

The  elements  of  the  mth  row  of  A  can  be  written  as 

(W{,{N~m  + '),  W{,(N-m+2),  ■  ■  •,  WkaN~m)),  (D4) 

or  using  Eq.  (D3)  equivalently  as 

(Wf,  (-m  +  I>,  wti~m+2)  ,  •  •  -,  (D5) 


If  am  is  a  vector  of  length  N  equal  to  the  mth  row,  then  the  condition  that  the  code  be  perfect  is 

a,aj  =  0,  i  ^  j ,  /,  j  =  I,  2,  •  ■  •  ,  /V,  (D6) 

where  t  is  the  conjugate  transpose  operation.  Note  that  Eq.  (D6)  is  representative  of  (N  -  IW/2 
equations.  However,  because  of  the  rows  are  rotations  of  one  another,  only  /V  -  I  equations  are  not 
redundant.  These  distinct  equations  can  be  written  as 

aja^  =  0,  m  #  1 .  (D7) 
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From  Eqs.  (D5)  and  (D7),  we  can  show  that  the  N  -  1  nonredundant  equations  are 

N-  1 

£  >-/<"-"»  =  0;  m  =  1,  2,  N  -  1.  (D8) 

n  =o 

We  know  that  if 

f  (n)  -  f  (n  —  m)  —  Mnm  +  g(m).  (D9) 

where  g  (m )  is  some  integer  function  of  n ,  and  M  is  an  arbitrary  integer  relatively  prime  to  N .  then 
Eq.  (D8)  is  satisfied  for  all  m  .  We  seek  the  functional  form  of  / ( n ). 

We  let  f(n )  be  a  polynomial  of  order  L : 

L 

/(«)  =  ,  (DIO) 

1=0 

where  the  c;  are  to  be  determined.  It  is  elementary  to  show  that  for  f(n )  having  the  form  given  by 
Eq.  (DIO)  that  in  order  to  satisfy  Eq.  (D9), 

ct  =  0  for  I  >  3.  (D1 1) 

Hence  f(n)  is  a  polynomial  whose  form  is  given  by 

f(n)  =  c,n  +  c2n2.  (D 1 2) 

(Note  that  we  have  dropped  the  c0  term  since  it  is  merely  an  arbitrary  linear  phase  shift  across  all  ele¬ 
ments.)  Now 

/ (n  —  m)  =  C|(n  —  m)  +  c2(n  -  m)2 

=  -2  c2nm  4-  cxn  +  c2n2  -  c  ,m  4-  c2m2 .  (D13) 

Hence,  if  we  subtract  Eq.  (D13)  from  Eq.  (D12)  and  use  the  condition  given  by  Eq.  (D9),  we  find 
that 

c2  =  j-M,  (D14) 

g(m)  =  ctm  -  c2ml. 

As  a  result, 

/(n)  =  Cjn  4-  Mn 2.  (D 1 5) 
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The  form  of  c ,  can  be  determined  from  constraint  congruence  (Eq.  (D3)).  After  substituting  the 
functional  form  of /(•)  as  given  by  Eq.  (D15)  into  Eq.  (D3)  and  simplifying,  we  find  that 

~MN2  +  c  iN  =  0  mod  N .  (D17) 

It  is  required  that  the  expression  on  the  left  of  this  congruence  be  an  integer  and  divisible  by  N .  It  is 
easy  to  show  that 

•  For  N  even,  c,  =  /,  arbitrary  integer, 

•  For  N  odd,  M  even,  c  t  =  / , 

•  For  N  odd,  M  odd,  c{  =  I  + 

Hence  we  can  write 


f(n)  =  ~Mn 2  + 


I  +  mod  2] 


(D18) 


If  we  drop  the  arbitrary  linear  phase  shift  across  all  elements,  it  can  be  shown  that  an  equivalent  form 
is  given  by  the  expression 

f(n)  =  yMn(n  +  N  mod  2).  (D19) 

The  P4  code  (8)  and  its  permutation  set  is  the  same  as  this  code,  which  we  call  a  P4/Chu  code, 
except  for  a  linear  phase  shift. 

REFERENCE 

Dl.  D.C.  Chu,  “Polyphase  with  Good  Periodic  Correlation  Properties,”  IEEE  Trans.  Inf.  Theory 
18,  531-532  (1972). 
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PROOF  OF  THE  ROTATIONAL  INVARIANCE  OF 
THE  APERIODIC  ACF  OF  THE  P4/CHU  CODE 


The  proof  for  the  rotational  invariance  of  the  P4  code  is  shown  here  in  terms  of  the  more  gen¬ 
eral  P4/Chu  code  described  in  Appendix  D.  The  general  form  of  the  P4/Chu  code  as  given  by  Eqs. 
(Dl)  and  (D19)  is 


an  =  wtf'Mnin  +N  mod  2) ,  „  =  0,  I ,  •  •  •  ,  N  -  1 


where 


.2  T 


WN 

M 

N  mod  2 


-j- 


1S 


an  integer  relatively  prime  to  N ,  i.e.,  (M ,  N)  =  1, 
fo  if  N  even, 

1  1  if  N  odd. 


An  arbitrary  rotation  of  this  code  that  forms  a  new  code  is  given  by 
a  (k)  _  +*)(n  +*  +N  mod  2) 


n  =  0,  •  •  • ,  N  -  1 
k  =  0,  1,  ,  N  -  1 


where  k  indexes  the  shifts  in  the  rotation.  We  define  a  code  word  a(t)  of  length  N  as 

The  aperiodic  ACF  sequence  is  defined  as  a  vector  r{k)  of  length  2N  -  1: 

r  -a  *a  -  — <vv  - 1}  •.  r-,w-2)>  >  ro  >ri  >  -rN-i)- 


(El) 


(E2) 


(3) 


(E4) 


where  *  denotes  linear  convolution  and  -  denotes  the  time  reversal  complex  conjugate  of  a.  We  can 
show  that 


N  -  1  -m 

.(k) 

m  L* 

n  =0 


Ak)„(k) 


m  =  0,  1,  2,  •••  ,  N  -  1 


(E5) 


and 


r(*)  _  _(*)* 
'  -m  ' m 


(E6) 


We  show  that  Ir^’l2  is  independent  of  k  or  equivalently  that  the  ACF  of  the  P4/Chu  code  is 
invariant  with  rotation. 
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To  this  end,  we  substitute  Eq.  (E2)  in  Eq.  (E5): 


.V  —  1  —  m 


r(k)  _  ^  +k){n  +k  +N  mod  2)  yy^-(\/2)M(n  +k  +  m)(n  +k  +m  +A/  mod  2)  (E7) 


n  =0 


New 


M(n  +  k)(n  +  k  +  N  mod  2)  -  M{n  +  k  +  m)(n  +  k  +  m  +  N  mod  2) 


=  —  2 Mkm  —  Mm 2  —  Mm  —  Mm  (N  mod  2)  —  2Mmn  . 


Thus 


N  —  1  —m 


=  w-vrnu&m+m'+m+mwmodl))  £  W„Mmn,  m  =  0,  1 ,  •  •  ■  ,  N  -  1 .  (E8) 


n  =0 


We  see  from  Eq.  (E8)  that 


N-l-m 

I  rm)  I  2  =  I  £  ^ 

n  =0 


(E9) 


Hence  the  magnitude  of  the  ACF  is  independent  of  k  or  invariant  with  rotations  of  the  original  code. 


i 

I 
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Appendix  F 

PROOF  OF  THEOREM  5 


Theorem  5:  Consider  an  N  x  N  matrix,  E  with  elements  on  the  unit  circle.  Let  E  have  mutually 
orthogonal  rows.  In  addition,  let  all  rotations  of  any  two  columns  of  E  be  mutually  orthogonal. 
Then  a  perfect  periodic  code  results  from  concatenating  the  rows  of  E. 


Proof :  We  define  cnm  to  be  a  vector  equal  to  the  «th  column  of  E  where  this  column  has  been  circu¬ 
larly  rotated  m  times.  To  be  a  perfect  periodic  code,  we  require  that  the  (N2xN2)  circulant  matrix 
T ,  whose  first  row  is  formed  by  concatenating  the  rows  of  E,  have  mutually  orthogonal  rows.  The 
matrix  T  is  given  by 


*11  *12 

’  *  1 N 

*21 

'  '  *iV  1 

eNN 

eNN  *  1 1  '  ' 

*  1  ,N  -  1 

*  \N  ’ 

eN  —  1  ,N  " 

eN,N- 

(FI) 


J_  *12  *  13  "■  *21  *22  eN  2  *1 


The  inner  product  between  rows  i  and  j  depends  only  on  i  -  j .  We  denote  the  inner  product 
between  row  1  and  row  m  +  1  by  rm  and  require  that  rm  be  equal  to  0  for  all  m  not  equal  to  zero. 

For  ms  0  mod  N ,  because  the  rims  are  mutually  orthogonal,  rm  =  0.  For  m  ^  0  mod  N ,  it 
is  straightforward  to  show  that 


o 

II 

+  (c:o.c10)  ' C  t, ,.C in )  +  .  .  . 

+  (C.voT.v-  |,o) 

-  (C|0,C^_! 

,|)  +  (C2„.Cvi,  +  (c  ,().C ,())  + 

■  •  ■  +  (c/VO'C/V-2.o) 

(F2) 


rV-’-1  ~  (c10'C2Ar)  +  (c20'C.w)  +  (Ci0.c4,V )  +  +  (CvOCl,/V-l)  < 

where  (•,•)  denotes  the  inner  product  of  two  vectors.  For  any  inner  product  term  (cn  ,  ckk  )  given 
in  Eq.  (F2),  n  j  ^  kx.  Hence,  if  all  rotations  of  any  two  columns  of  E  are  mutually  orthogonal,  it 
follows  by  inspection  of  Eq.  (F2)  that  rn  =  0.  Thus  the  theorem  follows. 
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Appendix  G 


PROOF  THAT  FRANK-P4  PRODUCT 
CODE  IS  A  PERFECT  PERIODIC  CODE 

The  Frank  code  of  length  N2  is  defined 

a(n)  =  Wfi {n/N],  n  =  0,  1,  N2  -  1  (Gl) 

where 

WN  =  e~Ji2T/N\  (G2) 


and  [  •  ]  is  the  least  integer  function. 

The  P4  code  of  length  N2,  where  N  is  odd,  is  defined 

b(n)  =  n$/2W*i(n  +  i)  t  n  =  o,  1,  •  •  ■  ,  N2  -  1  (G3) 


where  M  is  relatively  prime  to  N . 

We  form  a  new  code  called  the  Frank- P4  product  code  as  follows: 

c(n )  =  w$/2)Mn(m  +  l)  ■  W^n/N]  ,  n  =  0,  1,  ■  •  • ,  N2  —  1 .  (G4) 

We  can  show  that  c(n)  has  zero  sidelobes  when  implemented  as  a  periodic  code  if  M  +  1  is  also 
relatively  prime  to  N .  To  this  end,  by  using  Eq.  (G4)  it  can  be  shown  that  the  periodic  ACF  of  c(n) 
is  given  by 


N7-  1 

r(k)  —  ^2  +*)Kn  +*)//vl  jy(l/2)M(n  +*)(« 

n  =0 


+£  +  1) 


w~n[n/k\  yyO.Jl)Mn{n  +  1) 


(G5) 


k  =  1,  2,  ••  •,  N  -  1. 


We  can  show  that  Eq.  (G5)  reduces  to 


N7-  1 

r(k)  =  W<k7  *  ^2  +  G/N|  -  "[n/N]  yyMnk 

n  =0 


k  =  1,  2,  ■  •  •  .  N  -  1. 


(G6) 


Let 


n  =  n  ,/V  +  n2  , 

(G7) 

k  =  k\N  +  k2. 

(G8) 
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where  nx,  «2,  ku  k2  <  N .  We  rewrite  Eq.  (G6)  as 

r(k)  =  W(k*+k)M  ^n,'V+n-+*|N+Aj)I("|N+n!+*'W+*2)/yV1_<n,A/+":) l<nl/V+n:)/N] 

n2=0  n,  =0 

Simplifying  Eq.  (G9)  gives 


r(k )  =  W$(*2+<:) 


E  E  "V 

n,=0  n;=0 


u/Mnk 

"V 


Simplifying  further  yields 

r(k)  =  W$kl+k>W*N'k*  N£  *£ 

n,=0  n2=0 


-  <«■**>  <*■*■  21  e' 

n,=0  n,=0 


wMWtk)  "j,1  V' 


=  2'  EVr*""'*'  w?*""*'  <v- 


=  <*'*■  *£ 
n2= 0 

Now  for  k2  *  0,  and  A/  +  1  relatively  prime  to  N, 

*-•  (A 

E  ^ 

n.  =0 


A/-1 


E  ^ 

n,=0 


(Af  +  1)* 


:n,  1 


=  0. 


For  k2  =  0,  kx  #  0, 


r(k)  =  r(k,  1 V)  = 


e'  *  <w+u"2*' 

n2  —  0 


=  0. 


Thus 


r(k)  =  0  for  k  =  1 ,  2,  ■  ■  •  ,  N2  -  1 . 

Hence  the  Frank-P4  code  is  a  perfect  periodic  code. 


(G9) 


(G10) 


fG  11) 


(G12) 


(G13) 


(G14) 


Appendix  H 

PERMUTATION  CODES 


Consider  a  perfect  periodic  code  (PPC)  word  a  with  code  elements  an  n  = 
0,  1.  2,  •  •  •,  N  -  1.  Let  =  an.  Next,  consider  a  periodic  code  a'  with  code  elements  a'n, 

n  =0,  1,  •  •  •  ,  N  —  1,  where 


a 


n 


O Mn  mod  ,V 


(HI) 


and  M  is  relatively  prime  to  N .  We  show  that  a'  is  also  a  PPC  word.  Define  bm  such  that 


b 


m 


S-  1 


^  O  n  O  n  +  m  ■ 
n  =0 


(H2) 


Now,  a'  is  a  PPC  word  if  bm  =  0  for  m  -  1,  2,  •  •  •  ,  N  -  1.  We  know  that 

O  n  d n  *  m  +1  (H3) 

where  there  is  a  unique  k  and  /  for  each  n  and  m ,  respectively.  If  m  =£  0,  then  /  *  0.  Hence  a'  is 
perfect  since 


N-  1 

£ 


n  =0 


a'na 


'* 

n  +m 


N  -  1 


£ 

k  =0 


0*0?  +  / 


0. 


(H4) 


I 

1 
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INDEX  CODES 


Consider  the  following  code  with  P  elements: 

wk  '^d'<0, .  w[ in<v 1 ' ,  w[  ind■■'2, .  •  •  ■ ,  w[ ~ 1 1 


(II) 


where 


p 

is  a  prime  integer. 

a 

is  a  primitive  root  of  P , 

inda(- ) 

is  the  integer  index  function  to  the  base  a 

k 

is  a  divisor  of  P  -  1,  1  <k<P  -  1, 

l 

is  integer  not  divisible  by  n  , 

=  e  ~i2v/k ,  and 

w[ 

=  0  for  x  =0. 

We  show  that  the  periodic  code  given  by  Eq.  (II)  has  a  constant  sidelobe  level  -  1. 
We  can  show  that  the  mth  row  of  the  circulant  matrix  A  can  be  written  as 


I  ind  (P  +  1  -m  i  /  ind  (/>  +  I  ~m  +  I ) 

Wk  "  « 


WL 


I  ind  (P  *  I  m  ~  1 1 


(12) 


Let  bm  be  the  inner  product  between  the  first  row  and  the  m  +  1th  row.  then 


h-\ 

L 

n  =0 


w. 


I  ind  n  —  I  ind  In  -  ml 


k 


m  =  1. 


P  -  1 


(13) 


where  we  have  used  the  fact  that  inda  (P  —  m  +  n )  =  ind„(ti  -  m). 

A  theorem  [II]  in  number  theory  states  that  the  expression  given  by  Eq.  (13)  under  the  condi¬ 
tions  previously  given  for  P ,  a,  /  and  k  is  equal  to  —  1.  Hence 


bm  =  - 1,  m  =  \.  2,  •.  P  -  1 


(14) 


and  the  periodic  code  given  by  Eq.  (II)  has  a  constant  sidelobe  level  equal  to  -  1 . 
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